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MODELING THE RV AND BVS OF ACTIVE STARS 

Cezary Migaszewski ' and Grzegorz Nowak ' 



Abstract. We present a method of modeling the radial velocity (RV) measure- 
ments which can be useful in searching for planets hosted by chromospheri- 
cally active stars. We assume that the observed RV signal is induced by the 
reflex motion of a star as well as by distortions of spectral line profiles, mea- 
sured by the Bisector Velocity Span (BVS). The RVs are fitted with a common 
planetary model including RV correction term depending linearly on the BVS, 
which accounts for the stellar activity. The coefficient of correlation is an addi- 
tional free parameter of the RV model. That approach differs from correcting 
the RVs before or after fitting the "pure" planetary model. We test the method 
on simulated data derived for single-planet systems. The results are compared 
with the outcomes of algorithms found in the literature. 



1 Introduction 

The current spectral Doppler technique makes it possible to measure routinely the radial 
velocity (RV) of stars with a precision better than 10 ms^'. It is well known that the RV 
variability may be caused not only by the reflex motion of a star accompanied by smaller 
bodies, but also by internal activity of stellar atmosphere (e.g., non-radial pulsations, in- 
homogeneous convection, and/or rotating surface spots). While the reflex motion of the 
star leads to the Doppler shift of spectral lines, the atmospheric activity may cause distor- 
tions of spectral lines profiles (SLP), which displace their measured minima. Fortunately, 
a determination of the true origin of the RV variability is possible through the analysis of 
line profiles or the cross-coiTelation function (CCF). The basic tool for such analysis is 
the line bisector (LB) technique ( [Gray 1983]l2005l l. The most simple and useful measure 
of the slope of LB is the Bisector Velocity Span (BVS) defined as the difference between 
the LB velocity measured at some upper and lower flux levels of the SLR 

It is well known ( Desort et al. 20071 and references therein), that distortions of SLP 
caused by stellar activity may produce quasi-periodic RV signals mimicking planetary 
companions. The BVS analysis are then very helpful to reject or confirm the planetary 
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nature of the RV variability. The first case, when the BVS measurements were used 
to withdraw a hypothesis of a low-mass planet, was HD 166435 (Queloz et al. 2001) . 
In the next years, other stars with low-level RV variations stemming from stellar activ- 
ity were discovered, for instance, HD 78647 jSetiawan et al. 20041 . and HD 219542B 
jDesidera et al. 20041 . In the first two instances, the authors relied on a correlation be- 
tween the RV and the BVS data, indicating that periodicity of the RV signal has the 
stellar-activity origin. (Further in this paper we will demonstrate that, in general, such 
a correlation does not provide sufficient information to reject the planetary hypothesis). 
Recentlv. lBonfils et al. 2007l report a discovery of a "super-Earth" planet orbiting GJ 674 
with a 4.69-day period. Besides the strong signal of that planetary candidate, the RV 
data reveal also secondary, 34.85-day periodicity. The discovery team performed 2-planet 
Keplerian fits and, using photometric measurements, attributed the second period to the 
rotational modulation of the RV caused by a stellar spot. 

In this work, we propose an alternative method of resolving non-unique sources of 
the RV variability. It relies on simultaneous analysis of both sets of observables, RV and 
BVS, within an uniform fit model. In Section 2, we introduce the generalized model of 
RV. In Section 3, we simulate RV data for our experiments. Section 4 is devoted to RV 
models and the results, which are finally discussed in Conclusions. 

2 Modeling the RV and BVS data 

The stellar RV variability caused by the presence of additional bodies may be modeled as 
a superposition of Keplerian, astrocentric orbits jSmart 19491 ): 



where is the number of planets, Kj is the semi-amplitude of the contribution by the 
i-th planet, (53, is for the longitude of pericenter, e, is for the eccentricity, V; denotes the 
true anomaly (which depends implicitly on the orbital period P, and the time of pericenter 
passage T,), y„ are the RV offsets and O is their number. We search for such parameters 
{Ki,Pi,ei,(i5i,Xi,Vo), which may explain the RV variability in the sense of the least squares. 

When we deal with significant stellar activity, and the BVS measurements are avail- 
able, we can modify the RV model by adding a correction term, accounting for distortions 
of the SLP which contribute to the RV variabihty, = y,.(f) +AV{t). In the first 

approximation, AV{t) may be expressed through the BVS time series {{BVS} from here- 
after), AV{tj) — a{BVS}j, {j — 1, . . . ,N,n), where is the number of observations and 
a is a free parameter of the model. In fact, a depends not only on stellar activity but also 
on the spectral resolution (I Desort et al. 20071) . and a particular choice of the upper and 
lower segments of the LB. Now, following the principle of the least squares, we define 
the (Xr) '^■^ function of the fit model as follows: 




Vr{t) = £ Ki [cos (C3,- + V;) + e,- COS C3,] + Y^Vc 



(2.1) 



1=1 o=l 




(2.2) 
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Fig. 1. Simulated {RV} and {BVS} data for the following parameter tuples of 
(A'[ms"'],P[d],e,tD[''],3Wo['']); stellar-activity signal: (85,3.8,0.2,75,0); planetary signal: 
(150,3.8,0,270,180); stellar jitter is 30 ms"', Crv = 15 mi-^,<3Bvs = 10 ms"', a = -1.1364. 



where A^,, is the number of model parameters, Orvj, c?bvs.; stand for the standard errors 
of the {RV} and {BVS} time series, respectively. 

It is worth to note that the BVS may be used to correct the RVs for stellar activity 
contribution not only when the central star exhibits periodic variability, but also when 
that activity has an irregular (aperiodic) character. In fact, the BVSs measure shifts of 
detected minima of spectral lines, regardless of a type and sources of stellar activity. 



3 Simulated observations of the RV and BVS 



To give an example, and to demonstrate features of our method, we test model Eq. 2.2 on 
synthetic RV and BVS observations. We consider a star exhibiting periodic activity that 
mimics planetary RV signal, and we assume that it also hosts a planet. To construct the 
synthetic data set, we fix the elements of Keplerian orbit, and we simulate the planetary 



RV signal, {7?y}'^P'*, Eq. 2.1 We also simulate the {BVS} data. The synthetic signal, 
{RV} = {7?y }'P'^ + a{BVS}. Next, we add Gaussian errors independently to both sets of 
observables and also non-periodic noise (jitter) which is added to the {RV} and {BVS} in 
the same phase. Parameters of the synthetic signals are given in caption to Fig.[T] 

For a reference, we assume that a single-planet system has the same orbital and stellar- 
activity periods. Our choice of stellar parameters follows the first discovery of a star that 
mimics the planetary signature of RV, i.e., HD 166435 ( |Queloz et al. 200 1] ). The synthetic 
curves of the RV and BVS, as well as simulated measurements (the red points) in random 
epochs, are illustrated in Fig.[r| The first two panels are for the {RV} and {BVS} time 
series. The third panel is for a correlation between these observables. Grey, thick line 
illustrates the linear regression of the stellar-activity contribution to the RV {{RVY'^^^) on 
{BVS}. Black, thick fine marks the linear correlation between the {RV} and {BV^} sets. 



4 Tests of fitting algoritfims and the results 

The results of analysis of the synthetic data are illustrated in Fig. [2] Panels in this figure 
are for color-coded maps of (Xr)'''^ planes of selected fit parameters. Contours mark 
the standard confidence intervals of (la, 2a, 3a), respectively. The filled, crossed circles 
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Fig. 2. Maps of (x?)'''^ for the test case. Simulated data are presented in Fig. [I] The left column: 
the (Pft,ei)-plane; the right column: the (P/^j/fi) -plane. Each row is for a different fit algorithm 
tested in this paper (see the text for details). Filled, crossed circles mark the orbital parameters of 
the nominal system (given in caption to Fig.[T]l. 



mark positions of the nominal solution in the {P[,,e[,)- and {Pi,,Ki,)-planes (the left- and 
the right-hand columns, respectively). To compute the (Xr)^'^^-inaps, we proceed as fol- 
lows. We fix a particular point {x,y) in the selected parameter plane, and we search for 
remaining best-fit elements of the RV model. To search for the best-fit solution, we apply 
the Monte-Carlo method to choose initial conditions for the fast Levenberg-Marquard al- 
gorithm (see, e.g., Gozdzie wski et al. 20081) . Such fits are repeated for all points {x,y) of 
a discrete grid in the given plane of model parameters. 

Each row in Fig. |2]is for a different method of modeling the RVs. We start from 
the most simple (but generally wrong) method I that relies on fitting single-planet model 



(Eq. 2.1 1 to the {RV} signal, without making any use of additional information "hidden" 



in the BVS. The results are illustrated in Fig. [2^,b. Clearly, the nominal solution lies 
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beyond the 3a level of the best fit solution. In particular, the fitted Kh is significantly 
larger than the reference value. Both signals, {/?y}*P'^ and a {BVS}, have the same 
periods, and depending on their relative phase, the semi-amplitude of the resulting {RV} 
may be larger or smaller than that ones of the planetary signal. In the tested example, the 
{RV{ signal has larger amplitude than the {/?y}*P'*. Hence, we obtain too large mass of 
the planetary companion, however, the original orbital period is found correctly. 

To improve method I, one may proceed as follows (method II). At first, we find the 
coefficient ttobs of the linear correlation between the {RV{ and {BVS} signals. Then we 
subtract a correction term, aobslSV-S}, from the observed {RV} signal. We obtain, as 
we suppose, a pure planetary signal, so we can fit the single-planet model to corrected 
{RV} data. The results are presented in Fig.|2j;,d. Again, the derived best-fit solution is 
displaced from the true position (in particular, its semi-amplitude is badly determined). 
To explain that improper outcome of the fit, we look at the right-hand panel of Fig. [T] The 
grey line is for the correlation between {BVS} and {RV} observables when the star would 
be alone. If the planet is present then the detected correlation (ttobs, black line) changes. 
The difference between a and ttobs depends on the semi-amplitudes of both signals as 
well as on their relative phase. Unfortunately, also uncertainties of the measurements 
(instrumental ones or stemming from the stellar jitter), as well as their irregular sampling 
may change the correlation coefficient. Then, if we use aobs to correct the RVs by the 
linear term of (Xo\,s{BVS}, we may subtract too large or too small correction from the 
{RV} signal, and the best fit parameters of the orbital solution will be wrong. In contrary 
to method I, we obtain too small mass of the planet. 

The next approach relies on fitting 2-planet Keplerian model to the {RV} data (method 
III, see lBonfils et al. 20071) . We obtain two best-fit Keplerian orbits, moreover, we have 
to identify which orbit describes the true reflex motion of the star, and which one cor- 
responds to a false orbit mimicked by stellar activity. We did a test of this method and 
the results are illustrated in Fig. |2^,f. Still, we do not obtain correct results. The main 
difficulty in this approach emerges from a non-unique identification of the signals. For 
instance, a sum of two different (de-phased) RV signals with the same period may look 
very similar to a composition of two signals with the 2 : 1 commensurability of periods 
(see, e.g.. lGozdziewski and Konacki 2 006 ). In fact, the second orbital period found in the 
best fit solution, may be very different from its real value. Method III gives also too large 
mass of the planet. 

Finally, we apply our method IV relying on simultaneous analysis of both sets of 
{RV}- and {By5}-observables. We assume that a is an additional, free parameter of the 
fit model. Actually, we found that this approach provides the best results, as compared to 
the outcomes of methods I-III. That is illustrated in Fig.|2^,h. In this case, the agreement 
of the reference parameters with the orbital elements of the best-fit solution is basically 
perfect. Thanks to the self-consistent fitting process, we also find the true correlation 
coefficient a of the stellar contribution to the RV, {/?y }^"'', and the BVS. 

5 Conclusions 

Recently, the line bisectors are routinely derived from the same spectra used to measure 
the RVs (e.g., the Ii-cell technique, .Martinez Fiorenzano et al. 2005 J . They are much 
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more useful to understand the origin of the RV variations than alternative, even more 
indirect observables. The BVS data may help to detect stellar-induced periodicity in the 
RV signals or to correct the RVs by accounting for the linear correlation between the 
RV and BVS data. In this work, we also show that these data make it possible to find 
correct orbital parameters through fitting the planetary model to the RV and BVS data 
simultaneously. Our method seems to work well, even if there is a close commensurability 
between periods of both signals. In such a case, other algorithms tested here, provide the 
best-fit orbital elements which are significantly different from the true, reference solution. 

We note that the proposed method omits one more indicator of a distortion of the SLP, 
the Bisector Curvature (BC). The BCs measure the second order derivative of the LB. It 
still may be possible to improve the algorithm by using the BCs to correct AV{t). We 
are going to investigate such an improved method with more sophisticated and reahstic 
simulations of the measurements. 
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